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I. INTRODUCTION 
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A Quantum Kinetic Master Equation (QKME) for bosonic atoms is formulated. It is a quantum 
stochastic equation for the kinetics of a dilute quantum Bose gas, and describes the behavior and 
formation of Bose condensation. The key assumption in deriving the QKME is a Markov approx- 
imation for the atomic collision terms. In the present paper the basic structure of the theory is 
developed, and approximations are stated and justified to delineate the region of validity of the 
theory. Limiting cases of the QKME include the Quantum Boltzmann master equation and the 
\ Uehling-Uhlenbeck equation, as well as an equation analogous to the Gross-Pitaevskii equation. 
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£\J , Recent observations of Bose- Einstein condensation (BEC) 
in a dilute gas of magnetically trapped alkali atoms 
0-IHl-QSil' an< ^ m excitonic systems [?J has stimulated 
new theoretical efforts to describe the dynamics and sig- 
£2 ' natures of weakly interacting Bose gases. In contrast 
to superfiuid Helium, which for many years has been 
the only experimental example of BEC, these weakly 
y—i , interacting Bose gases are much more amenable to a 
^sO • theoretical analysis. Recent theoretical work has fo- 
0^ \ cused on a description of BEC in a trapping poten- 
^ tial and dynamics of formation of a Bose condensate 

In this paper we want to develop techniques to de- 
scribe the behavior and formation of the Bose condensate 
in which we will combine the simplicity of a quantum 
j—* ( stochastic methods used in quantum optics [^0[ with a 
realistic treatment of the interatomic interactions which 
are known to play a major role in the dynamics of the con- 
densing system pi) . Such a description must necessarily 
involve aspects of a kinetic theory as well as quantum 
C$ ' mechanics. 

The first realistic formulation of kinetic theory was the 
Boltzmann Equation (BE) whose usefulness and success 
to this day is universally accepted. The so-called Quan- 
tum Boltzmann or Uehling-Uhlenbeck equation (QBE) 
J42| introduces corrections for quantum statistics into the 
Boltzmann collision term, but since it deals only with 
one-particle distribution functions we cannot expect this 
equation to give a realistic treatment of the quantum me- 
chanical aspects which must occur in a Bose condensate. 

A different description is provided by the time depen- 
dent Gross-Pitaevskii equation (GP equation) |Q which 
can be viewed as an equation for the condensate wave 
function (order parameter for the Bose condensate) . This 
equation is clearly a simplified description in that it in- 
cludes no quantum fluctuations, or thermal or irreversible 
effects, but it may well be valid in the situation of a large 
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number of condensate particles. Both of these equations 
contain essential aspects of the problem we are study- 
ing. However, in practice the process of creating a Bose 
condensate in a trap by means of evaporative cooling 
starts in a regime covered by a kinetic equation and fin- 
ishes in a regime where the GP is thought to be valid. 
What is needed is unified description which covers the 
whole range — a true quantum kinetic theory. In this 
paper we will develop a quantum kinetic master equation 
(QKME) which is not a mere modification of the BE but 
extends its philosophy into quantum mechanics and de- 
rives a quantum stochastic equation of the kind similar to 
those derived in quantum optics. 

The QKME which we derive is thus a consistent ex- 
tension of the philosophy of the BE to a regime in which 
quantum mechanical effects are not a minor correction 
to classical results, but play a major role in the full de- 
scription of the system. 

In applying the QKME to the problem of Bose conden- 
sation, although quantum mechanics plays major role, 
the significant quantum aspects are restricted to a few 
modes, the remaining modes being able to be describe 
in a manner similar to the classical Boltzmann equa- 
tion. This is a situation which is very familiar in quan- 
tum optics, where we very often find that it is neces- 
sary to describe a few optical modes fully quantum me- 
chanically, while treating the rest of the systems as a 
heat bath. However, the methods necessary for quan- 
tum kinetic theory contain considerably more technical 
difficulty, largely because of the interaction between the 
atoms. The nearest optical analogy is multimode prop- 
agation in a medium with a strong nonlinear refractive 
index @. 

II. THE NEED FOR A QUANTUM KINETIC 
THEORY 

Our goal is to derive a quantum kinetic master equa- 
tion for a gas of TV weakly interacting Bose atoms. The 
master equation is formulated as an evolution equation 
for the iV-atom density matrix, giving a quantum me- 
chanical generalization of the Quantum Boltzmann equa- 
tion. This provides a fully quantum mechanical descrip- 
tion of the kinetics of a Bose gas, including the regime 
of Bose condensation. In particular, such an equation is 
capable of describing the formation of the Bose conden- 
sate. The quantum mechanical processes which must be 
included in this treatment are the atomic motion (trans- 
port) and coherent and incoherent interactions (colli- 
sions) between the atoms. In this section we give a qual- 
itative overview over the main results to be derived in 
later sections. 
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A. Essential elements of the theory 

There are three principal elements, as follows. 

1. Coherent and incoherent dynamics 

The derivation of a quantum master equation implies 
that the system dynamics can divided into into coherent 
and incoherent physical processes, where the incoherent 
dynamics is modeled as a quantum stochastic Markovian 
process. This builds on ideas, which have been developed 
and applied successfully in quantum optics as a method 
of describing coherent and incoherent processes simulta- 
neously (e.g. in laser theory). 

In a weakly interacting Bose gas, we will treat the col- 
lisions between the atoms resulting in a large momentum 
transfer as random and weak incoherent processes, re- 
sponsible for the "noise" in the system. On the other 
hand, forward or near forward scattering is, if the wave- 
length is sufficiently long, largely a coherent process, giv- 
ing rise to dispersive effects. 

2. Dynamics of the condensate 

The (coherent) dynamics of the Bose condensate is 
treated explicitly and separate from the description of 
the non-condensed modes which play the role of noise 
and feeding terms for the Bose condensate. Again, this 
is guided by laser theory where the coherent laser mode is 
separated from the weakly populated incoherent modes 
which are typically treated as a heat bath. 

3. The "cell" description 

The theory is formulated in terms of a quantum me- 
chanical phase space description for with coarse grained 
position and momentum variables, based on division of 
phase space into "cells" of volume h 3 . This seems to 
be the most natural formulation, since collisions are best 
described in terms of momentum, whereas transport is 
best understood in terms of position. We will introduce 
this formulation by means of a wavelet description, which 
gives an exact exact description of the dynamics; approx- 
imations appear as a result of the procedure for describ- 
ing the kinetics, and the conditions for the validity of our 
approximations then put some conditions on the precise 
nature of the cells chosen. 



B. The scope of this paper 

In developing the physical picture, there are three main 
stages through which we must proceed in order to achieve 
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a full quantum kinetic theory applicable to the present 
experiments on Bose-Einstein condensation. 



i) It is necessary first to develop the theory for the weakly 
interacting Bose gas with no trapping potential, whose 
stationary solution is well approximated by an ideal Bose 
gas. 

ii) The next stage is to take account of the modifications 
(these may be very significant) to the dynamics, which 
are induced by the presence of a large proportion of atoms 
being in the condensate. This can be developed in terms 



of Bogoliubov quasiparticles jl5 46 1 , or more generally. 

iii) Finally, we introduce the trapping potential. There 
are two extreme situations to consider. If the trap is 
very broad, then the situation is qualitatively not very 
different from the case of no trap, though the appropri- 
ate modifications can be quite intricate. In the case of 
a very tight trap, we have a completely different treat- 
ment, based on the kinetics of transitions between differ- 
ent trap levels. There will also be an intermediate regime, 
in which it is relevant to treat only the lower trap levels 
explicitly, while upper levels are thermalized. 

This paper deals with i) only; ii) and iii) will be treated in 
our second paper QKII. Thus this paper develops all the 
basic ideas and methodology, while QKII will develop the 
refinements necessary to deal with strong condensation 
with and without a trapping potential. 



1. The physical picture 

This phase space description is one in which particles 
are represented as wave packets interacting with each 
other. From the point of view of wavepackets there are 
the following processes: Transport: that is the motion 
through space of a wavepacket. However we also under- 
stand that from a wavepacket point of view, transport is a 
process by which the wavepacket moves with unchanged 
shape — the wavepacket spreading inherent in the quan- 
tum mechanics of particle motion is not itself viewed as 
a part of transport. Wavepacket spreading occurs in ad- 
dition to what one would normally call transport, and 
must contribute to the production of coherence which 
can arise in a Bose condensed system. It is purely quan- 
tum mechanical, but arises of course from exactly the 
same quantum mechanical source as transport, as de- 
fined here. Motion in a trap is readily included in these 
arguments to the extent the trapping potential can be 
considered as slowly varying. Collisions: here we mean 
localized events in which the momentum of wavepackets 
will change instantaneously, with a consequent change 
in the wavefunction. The localization cannot of course 
be exact, because wavepackets themselves are extended 
objects. 
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2. The size of the cells 



Quantum mechanics itself does not give any preferred 
size of the momentum space or coordinate space cells 
separately. Physically, we can characterize the choice 
of cells by a wavenumber A (or a cell size l c = n/A 
(compare Fig. 1), and this will naturally provide a "cut" 
in momentum space; changes of momentum greater than 
TiA are seen as collisions, while smaller changes are seen 
as changes which take place in coordinate space. The 
value of A is constrained by two considerations. 

(a) We want to have a description in terms of "en- 
ergy levels" within a box inside each cell. This 
means that the wavefunction of a particular 
particle must maintain its coherence over the 
length l c of the cell. The coherence length 
which arises from collisions should normally 
be of the same order of magnitude as the mean 
free path A m f p = vt, where v is the speed of 
the particle, and r is the mean time between 
collisions. Thus this requirement yields the 
condition l c <C A m f p . 

(b) We will be wanting to develop a Markovian 
description of collisions, which depends on the 
determination of a characteristic decoherence 
time, which is essentially the bandwidth of fre- 
quencies of particles in the system, that is, 
thermal correlation time h/2irkT. This should 
be much less than the the characteristic time 
for evolution of the the phase-space density, 
and this yields the condition that the mean 
free path A m f p 3> At = h/y/2mkT, which is 
the wavelength of a typical particle, known as 
the thermal wavelength. 

(c) As well as the condition that the bandwidth of 
frequencies must be sufficiently large, we must 
also have a sufficient density of states to ensure 
that the correlation function inherent does in- 
deed become very small for times greater than 
the correlation time. Essentially this requires 
that the maximum energy difference between 
cells should be much less than kT, lea ding to 
the requirement l c 3> It = h/V2mkT, that is 
the cells much be very much larger than the 
thermal wavelength. 

(d) The method is based on a requirement that 
correlations which involve momentum differ- 
ences greater than A are negligible. This 
means that the "collisions", that is, processes 
which cause a momentum change larger than 
|A|, are seen as inducing incoherence. The 
derivation of the QKME which we use does 
this by means of a projection formalism, and 
the condition for the validity of the procedure 
is found to be 
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aA m fpAT <C l c 



(1) 



where a is the scattering length for the inter- 
action between the particles, 
(e) There is also a weak condensation condition 
for the validity of the methods used in this 
paper, and this can be written in terms of the 
cell size l c , the scattering length a, and the 
particle density p as l c ^> yjir /8ap. 

These criteria are easily met, for example in the ex- 
periment of H one finds that A = 0.4619m, while the 
thermal wavelength is It — 4.8 x 10~ 7 m. 

On the other hand in that experiment the trap itself 
has a size of the order of magnitude as It- Thus to de- 
scribe the behavior of the wavefunction within the trap, 
one would want l c to be finer than this. This first paper 
does not deal with trapping potentials, but we will show 
how to overcome this difficulty in the QKII. The essence 
of the solution is that it is not necessary to have the same 
size cells for all ranges of momentum — they can be finer 
in space at the lower momentum ranges. 

Figs. 2 and 3 illustrate the phase space description 
and the various physical processes. The figures show the 
phase space cells corresponding to different momentum 
bands K with width A and position cells r. A state 
of the system is given by specifying a set of occupation 
numbers {n(K,r)} of the cells (Fig. 2). Transport and 
wave packet spreading corresponds to coherent quantum 
processes connecting cells in a given momentum band 
(horizontal arrows in Fig. 3a). Note that these coherent 
processes do not change a given momentum distribution 
N(K) = (-/V(Ki), 7V(K 2 ), . . .) with N(K) = £ r n(K,r) 
the number of atoms in a momentum band. In a sim- 
ilar way, coherent processes corresponding to a forward 
scattering, or motion in a slowly varying trapping po- 
tential give rise to smooth shifts between the cells of 
neighboring K (Fig. 3a and b). Incoherent collision pro- 
cesses between two atoms, on the other hand, are asso- 
ciated with a "quantum jump like" momentum change 
Ki,K2 — > K3,K4 between widely separated in momen- 
tum bands (however, essentially confined to within a sin- 
gle spatial cell r) (see Fig. 3c). These transitions conserve 
momentum, 

K! + K 2 = K 3 + K 4 + 0(A), (2) 

up to an uncertainty A associated with the width of the 
band, and energy 

E 1 +E 2 = E 3 + E 4 (Ei = h 2 K 2 /2m). (3) 

Collisions cause transitions between different momentum 
distributions 

N(K) = (. . . ,71l, 71 2 , 71 3 , 71 4 , ■ ■ ■) 

-» N'(K) = (. . . , m - 1, n 2 - 1, n 3 + 1, n 4 + 1, . . .) (4) 
which we write in the form 

N(K) — > N'(K) = N(K) + e. (5) 
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C. The quantum kinetic master equation 



The derivation of the master equation in Sec. [V C be- 



low assumes that the full quantum coherence is kept for 
different spatial configurations (occupation numbers of 
the phase space cells) within one K-band, and the as- 
sumption that atomic coherences between different mo- 
mentum bands are eliminated in the Born and Markov 
approximation. The mathematical procedure of deriv- 
ing a master equation in Sec. 
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is based on projecting 
the iV-atom density operator p on states associated with 
a given momentum distributions N(K). For a given cell 
size there is a regime where this diagonally projected den- 
sity matrix wn obeys a closed equation 



*n(*) = -- 



E 

. K 



I 

2h 



n 2 v 2 

2m 



V'k(x), Vys(t) 



(C/(K 1 ,K 2 ,K 1 ,K 2 ) + f/(K 1) K 2 ,K 2 ,K 1 ))+^f/(K 1) K 1 ,K 1 ,K 1 ), UN (t) 



Ki#K 2 



-^£(A£(e)){2£/(e)« N _ e (t)£/t(e) - uHe)U(e)v N (t) - v N (t)tf(e)U(e)}. 



(6a) 

(6b) 
(6c) 



The first line in (f3aj) describes atomic motion between 
the different phase space cells within a given K-band. 
The first term is a streaming term with current operator 
j K (x) (for a formal definition see (^ below), and the sec- 
ond term corresponds to wavepacket spreading (the field 
operator for a K-band will be defined in Eq. (eql.16)). 
The second line is (coherent) forward scattering, and the 
last describes the redistribution by collisions according 
to a [/(e)-operator, which is expressible in terms of the 
interaction potential, and is explicitly defined in ([l7|). A 
detailed derivation of the equation will follow in Sec. IV 
below. 



1. Wavefunction stochastic differential equation (SDE) 
interpretation 

It is interesting to interpret ( |6a[ ) from the point of view 
of evolution of a stochastic iV-atom wavefunction where 
the wide angle collisions are described as quantum jumps 
|{47f . Note that ( |ua| ) can be written in the form 

»**(*) = -i^effWN(i) + lv N {t)Hl s (7) 

+ T E S ( AE ( e )) U(e)v N - e (t)U^e) . (8) 

e 

Let us define an atomic state vector \ip^(t)) with N(K) 
a given configuration of occupation numbers. In the time 
evolution of this state vector a collision is associated with 
an instantaneous jump 

\ip N , +e (t + dt)) cx U{e)\ips(t)) (9) 
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Here the [/-matrix plays the role of the jump operator, 
connecting configurations N(K) — ► N(K + e). The prob- 
ability for a collision is in the time interval (i, t + dt] is, 
given a normalized statevector IV'n(k)W) a t time t, by 

*W] = T ^5(Ai?(e))||[/(e)|^NW)|| 2 d« (10) 

e 

= ^P ( t, t+dt ](e) (11) 

e 

The time evolution between the collisions is governed by 
the nonhermitian Hamiltonian H e s, 

|^N(K)(t))=e-" r - tt / R |^(K)(0)). (12) 

Physically, (jL2h describes the streaming and forward 
scattering between the jumps, and the nonhermitian is 
associated with loss due to collisions. Note that H e g 
preserves the distribution N(K) (see Figs. 2a and b). A 
stochastic average over these quantum trajectories give 
the density matrix, v N (t) = ((|^ N (t))(^ N (t)|/||^N (t)| 2 )). 



D. Limiting cases 



The Quantum Kinetic Master Equation (6a) includes 
as a limiting cases both the Quantum Boltzmann master 
equation and the Uehling-Uhlenbeck equation, as well as 
an equation analogous to the Gross-Pitaevskii equation. 
These results will be derived in Sec. |v|). 



1. Quantum Boltzmann master equation 

If we assume a single spatial cell, i.e. the size of the 
system equals the dimension of the system (compare 
Figs. 1 and 2), the transport terms trivially disappear 
and Eq. (6a) reduces (Sec. [v]) to an equation which we 



shall call the Quantum Boltzmann master equation: 
u>n = ~^J2 <5(A£(1234))|Ti 234 | 2 

1234 



x|nin 2 (n 3 + l)(n 4 + l)[w n - w n +e] 

+ (ni + l)(n 2 + l)ri3«4[w n - Wn-e]}- (13) 



This is a rate equation connecting the occupation prob- 
abilities of a given configuration n = N, u>„, where the 
factors 7ifc + 1 reflect quantum statistics for the collisional 
transition rates. It obviously has a limited validity, but 
it is attractive because of its simplicity and ease of sim- 
ulation. El. 



2. Uehling-Uhlenbeck equation 
For spatially inhomogeneous systems w e de rive from 



(|6a|) the Uehling-Uhlenbeck equation (Sec. ( VB)). Let us 
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define a single particle distribution function /k (x) where 
K labels a momentum band and x is a (contin uous ) spa- 



tial coordinate (for a precise definition see Eq. ( pL 1 2[ ) ) . We 
interpret /k(x) as a joint momentum position distribu- 
tion function, similar to the Wigner function. Assuming 
a factorization of the TV-atom distributions (an assump- 
tion not valid in the BEC regime) one obtains from j6a| ) 
the kinetic equation 

| /K (x)) » ^/kW + ****** 

x6(K + K 2 - K.3 - K 4 )5(cj + w 2 - ^3 - W4) 
x{/k(x)/k 2 (x)[/k 3 (x) + 1][/k 4 (x) + 1] 
-[/k(x) + 1][/k 3 (x) + 1]/k 3 (x)/ K4 (x)} 

(14) 

which now includes a streaming term. The collisional 
term in Eq. ( jl4| ) as been approximated as s-wave scat- 
tering. 



3. Condensate master equation 



Finally, to illustrate t he k inetics of Bose condensation 
we will derive in Sec. ( |VCD a simple approximation of 
the QKME based on the assumption that the density 
operator of the total system can be factorized into a con- 
densate density operator po(t) for the K = Bose con- 
densate band, and an operator for the non-condensed 
modes K / which is assumed to be in thermal equi- 
librium. We obtain the equation, which we shall call the 
condensate master equation: 



Po(t) 



drxipl (x 



ft 2 V 2 
2m 



ipo (x) 



1 



rf 3 xV'( i )(x)^( i ,(x)-0o(x)'0o(x) , po 



(15) 



( Uff (0)^n K ) / d 3 x^(x)^ (x), 



d 3 x / d 3 x'G(-)(x-x',T, At ){2^ (x)p ^(x') -poV'i(x / )V'o(x) - Vj(x')^o(x)p } 



d 3 x J A'G(+»(x- x',T,/i){2V4(x)poV>o(x') - p ^o(x')^(x) - Vo(x')^(x)p }, 



(16) 



Here V'o(x) refers to the atomic destruction operator for 
the K = band with the Bose condensate. The first line 
in Eq. (|l^) gives the dynamic of the Bose condensate 
including the nonlinear interaction term proportional to 
ipQipQipoipo > an d an interaction term of the atom in the 
Bose condensate with the above condensate particles 
with thermal occupation numbers g(0) "-k^o (where 
<?(0) is a normalization factor to be defined in Eq. (p9[)). 
For zero temperature, and when the field operators "00 ( x ) 
are replaced by c-numbers, Eq. (15) is, of course, equiva- 
lent to the Gross-Pitaevskii equation. For finite temper- 
ature, solutions of the corresponding Hartree Fock equa- 
tions have been discussed in references [^9[|0 21 . The 
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second and third line in Eq. ( |l5| ) give dissipative loss and 
feeding terms for the condensate due collisions with the 
non-condensed atoms. The relevant collisions are of the 
form Ki, K2 •*-> K3, K = 0, where atoms are transferred 
to the Bose mode under conservation of momentum and 
energy. In Eq. (Ill]) these transition rates are denoted to 



GW(x-x', T, h) (to be defined in Sec. |3§); they involve 
the occupation numbers of the K ^ modes according 
to the given temperature and chemical potential 



4- Validity of the limiting cases 

These three cases are all very simplified, and more re- 
alistic applications of the QKME will require much more 
care. But each one them illustrates a different aspect of 
the problem of quantum kinetics, which gives them great 
value as an aid to intuition. 



III. DESCRIPTION OF THE SYSTEM 

We consider a set of Spin-0 Bose particles, described 
by the Hamiltonian 



H = Hq + Hi + H T , (17) 



in which 



V>(x) (18) 



and 

Hi = ^Jd 3 xJ d 3 xV t (x)V' t (x')w(x-x')V'(x')V'(x). (19) 

Thus, the operators ^>(x) have the commutators 

h&(x), VV)]=(S(x-x'), (20) 

[V(x)^(x')] = [V t (x),^(x')] = 0. (21) 

The potential function u(x — x') is a c-number and is 
as usual not the true interatomic potential, but rather a 
short range potential — approximately a delta function — 
which reproduces the correct scattering length. This en- 
ables the Born approximation to be applied, and thus 
simplifies the mathematics considerably. |50| 

The term Ht arises from a trapping potential, and is 
written as 

H T = J d 3 xV T (x)^(x)^{x). (22) 

In this paper we will restrict our considerations to the 
case of no trapping potential, so Vr will be set equal to 
zero. The modifications necessary if there is a nonzero 
trapping potential will be dealt with in a future publica- 
tion. 

Thus, this is the standard second quantized theory of 
an interacting Bose gas of particles with mass m. 
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A. Possible phase space descriptions 



1 . The Wigner function 

The original work of Bose and Einstein was con- 
ceived before the invention of quantum mechanics by 
Heisenberg and Schrodinger, and was based on a simple 
phase space description, in which both momentum and 
coordinate space were discretized to give phase space cells 
of volume h 3 . The state of the system was then speci- 
fied by the number of quanta in each phase space cell. 
The treatments in modern elementary textbooks use es- 
sentially the same arguments, but use a single large box 
in coordinate space, so that the discretizing comes about 
from the discrete energy levels in this box. 

We will want to develop a formulation which is as close 
to the Boltzmann equation as is permitted by the quan- 
tum mechanical nature of the problem. This requires 
a phase-space description of the field ?/>(x). Although it 
might seem that the Wigner function provides the appro- 
priate method, it is in fact impossible to write a phase- 
space Wigner function for a multiparticle system which 
manifestly exhibits the symmetry of the Bose wavefunc- 
tion. A simple demonstration of this problem can be 
made for a two particle system whose Wigner function is 
obtained from the Fourier transforms with respect to yi , 

V2, of 

Tr{V>(xi + yi)V>(x2 + y2)^ f (xi - yi)V' t (x2 - y 2 ) i o } 

= W(x 1) yi,x a ,y 2 ), (23) 

where po is the density operator of the vacuum. Thus, 
the Wigner function is (where TV is an appropriate nor- 
malization factor) 

W r (x 1 ,k 1 ,x 2 ,k 2 ) = 

N J e 2 ^ kl - yi+k2 - y2 )VK(x 1 , yi ,x 2 , y2 )d 3 yid 3 y2. (24) 

The Wigner function has the symmetry 1 <-» 2 ; but this is 
not the full Bose symmetry, which requires full symmetry 
under the independent exchanges: 

xi+yi/2^x 2 +y 2 /2, (25) 
xi - yi/2 <-> x 2 - y 2 /2 . 

The Wigner function symmetry 1 <-> 2 is given by the 
simultaneous exchanges (p5|), not the independent ex- 
changes. The problem obviously persists for all numbers 
of particles. Where the Wigner function has been used in 
kinetic theory J52] the Bose symmetry takes the form of a 
very complicated integral operator relation, whose use is 
very impractical. Although an exact evolution equation 
would preserve Bose symmetry, the absence of manifest 
Bose symmetry makes it extremely difficult to develop 
approximation methods which preserve Bose symmetry. 
In fact it is only in a very dilute gas limit, in which essen- 
tially only the one-particle Wigner function occurs that 
the Wigner function has been used, and for this of course 
Bose symmetry is not an issue. 
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2. Wavelet expansion of field operators 



The most natural way to give a phase space descrip- 
tion of a Bose gas is to follow the method of Bose and 
Einstein — that is, to divide space into cells of volume 
AV, and for each cell to introduce a set of basis wave- 
functions which vanish outside the cell, but are quan- 
tized within the cell. Unfortunately, even from the point 
of view of the non interacting Hamiltonian (|l8f), such 
wavefunctions have infinite energy, arising from the sharp 
transition from inside to outside the cell. This yields a 
wavcfunction which is not twice differentiable, as required 
by (|l^), so that the spread in momentum p is so large 
that the mean of p 2 is diverges. 

To avoid this problem, the transition from inside the 
cell to outside it must be made smoother, and this is a 
problem which is like that which arises in the study of 
wavelets |53| . We introduce a set of wavelet functions (in 
one dimension) 



f K+A 

v K {x,r) = . / e^-^dk (26) 

V4ttA Jk-a 



1 



(27) 

If 

r = mr/A, n = 0, ±1, ±2, ... (28) 
these functions have the property of orthogonality: 

dxV^(x,r)vK'(x,r) = SKK'^rr' (29) 



and they are also complete. From the definition (26), 
it can be seen that vx{x,r) has momenta in the range 
[h{K - A), %{K + A)), and from of @, it can be seen 
that wavefunction vk(x, r) is localized to a certain extent 
at the point x — r = mr/ A. 

These wavefunctions correspond to a complete set of 
phase cells; we see that if the uncertainties are defined as 
the intervals between the discrete values of the momenta 
and position, then 

5x = tt/ A and Sp = 2 Ah => Sx 5p = h (30) 

However, the uncertainty as defined by a variance in x is 
clearly infinite, reflecting the fact that the wavelet func- 
tions (|2^),(p7|) are not well localized. The quantities 8x, 
8p represent rather the spacing between the phase cells — 
not the uncertainties in x and p. In terms of the variable 
x, these wavelet functions represent the smoothest pos- 
sible functions we could choose, because by construction 
the momentum spread is bounded, so that the wavelet 
functions are in fact infinitely differentiable. 

The field operators can now be expanded as (in an 
obvious three-dimensional generalization) 
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^ x )-EE u K(x,r)a(r,K) (31) 

K r 

and the commutation relation for a(r, K) is 

[a(r,K),at(r',K')]=<W5 K K'- (32) 

The states of the Bose gas are now specified by the eigen- 
values of the number operators N(r, K), which is a truly 
quantum mechanical version of the original idea of Bosc. 

Notice that there is no need in this description for A 
to be unique — one can choose a different value of A for 
each K-band if one wishes, since momentum bands are 
orthogonal independently of their size. 



3. Momentum resolved field operators 



For many parts of our discussion, it is convenient to 
resolve field operators into only the different ranges of 
momentum; thus we write 



K 



where 



V.K(x) = e lK - x ]T^(x,r)a(r,K) 

r 

j L> K (x-x')V(xW 



iK-x 



where 



£> K (x - x') = e- 4K - (x - x,) .g(x - x') 



where 



7(x) = 



sin Aa 



sin Ay 



sin Az 



We also have the commutation relation 

[^k(x), Vk'( x ')] =<*kk'0(x-x'). 



(33) 

(34) 
(35) 

(36) 

(37) 

(38) 



The resolution into momentum bands thus yields a rather 
nonlocal description. Notice also that we can write 



<?(x-x') = 



1 



(27T) 



A 



s ik.(x-x') d 3 k 



^ e - K -( x - x V(x,r)^(x',r) 



(39) 



which is an expression of the completeness of the wavelet 
functions within a i^-band, and is useful in some of the 
computations u sed in deriving the Uehling-Uhlenbeck 



equation in sect. VB 
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4- Free Hamiltonian in terms of momentum resolved field 
operators 



Because the functions used in the expansion of e * K x V>k (x) 
are orthogonal to those used for e~* K x ^k'(x), it is ob- 
vious that 

yd 3 x^(x)^ K '(x)e- 4(K '- K) - x = ifK^K'. (40) 
Using this fact, it is straightforward to show that 

H ° = E / ^k(x)V>k(x) (41a) 

+ Y, /rf 3 x^K-j K (x) (41b) 

K 

+ E/ rf3x ^(x) (-^) V-k(x) (41c) 

= -ff a + fffc + i? c 

Here we have used a "probability current" for the mo- 
mentum K defined in a way analogous to that normally 
employed 

Jk(x) = ~ {<(r)V^ K (r) - V^(r)^(r)} . (42) 

The resolution into three parts has a simple interpreta- 
tion. The resolution into a full wavelet description shows 
that we can write 

^ = E^ a K««K(r) (43) 

corresponding to an energy ?i 2 K 2 /2m for each quantum 
at each phase space location K,r. This is the Hamil- 
tonian corresponding to the original ideas of Bose and 
Einstein. 

The term corresponds to transport, as shown by 
computing the commutator 

[H b , V>k(x)] = ikK • V^k(x) (44) 

where vk = KK./m. This is obviously transport corre- 
sponding to the velocity appropriate to the central mo- 
mentum of the momentum band K. 

Finally, the last term H c corresponds to wavepacket 
spreading through 

[Hc^k(x)] = -£v 2 Vk(x) (45) 

Thus the wavelet description naturally gives the separa- 
tion into the three parts fundamental to a phase space 
description of the processes. 
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5. Interaction Hamiltonian in terms of the momentum 
resolved field operators 



It is trivial that the interaction part, Hi, can be writ- 
ten 

Hi = \ J2 ^(Ki,K 2 ,K 3 ,K 4 ), (46) 



2 

Ki,K 2 ,K 3 ,K 



in which 



U(Ki,K 2l K 3l K4) = 

d 3 x I <i 3 x' g( iK l' x + iK 2- x '- iK 3-X-lK 4 -x') 



x A, WV£ 2 (x>(x - x')^k 3 (x')V-k 4 (x). (47) 
Furthermore, it is easy to see that 

[/(K 1 ,K 2 ,K 3 ,K 4 ) = (48) 

unless 

Ki + K 2 = K 3 + K 4 + 0(A) (49) 

i.e., the quasimomentum Ki is conserved to within some- 
thing not much bigger than A. 

IV. TREATMENT OF COLLISIONS 

A. Eigenstates of H a and its corresponding 
Liouvillian 

Suppose we use the wavelet ba sis; t hen we can label the 



eigenstates of H a (as given in (41a)) as |n) where n is 
a vector of elements rii , the number of quanta belonging 
to the phase space wavelet i, which has quantum number 
(Kj, Tj). The eigenstates of the corresponding Liouvillian 
are given by the outer product |n) (m| , and the eigenvalue 
is only zero if the energy of |n) is the same as that of |m) . 



B. Projectors 

We will define a projector onto all eigenstates with the 
same number of quanta in each K, irrespective of r; thus 
we define 

p N |n) = |n) if ^n(K,r) = N(K) and K ^ 0; 

r 

= otherwise. (50) 

This kind of projector identifies all configurations with 
the same distribution in K, but ignores K = and leaves 
undisturbed the distribution over r. 

Suppose now that p is the density operator for the 
system of particles; we define a projector on p by 
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V N p = pnPPn = un- (51) 

Thus the projected p is an operator in position space, 
but acts only with a space that has the configuration N 
of particles over the momentum space cells — except any 
configuration of K = particles is permitted (we call this 
K diagonal). It follows also that, for any p 

V^H bP = H b V N p (52a) 
VnH cP = H c VkP (52b) 

and that for H a there is the stronger result 

[H ai V N p]=0 (53) 

which is true by construction. 

The thrust of this paper is now to develop a closed 
equation of motion for the projected parts of p(t), namely 
Vjsi(t), and to show that we can use these parts to repre- 
sent the physics in which we are interested. Thus we will 
be making the approximation 

p(t)«£uN(*)> (54) 

N 

which will require the assumption that we can neglect 
the remainder: thus 



P(t) 

w 0. (55) 

A density operator of the form (|54|) will be called K- 
diagonal. For such if-diagonal density operators, the al- 
gebra in deriving the equations of motion is considerably 
simplified by results like 

Pm{4W^(x')%} = 

unless K = K' and M = N. (56) 

Notice however that in this equation, we do not require 
that x = x'; thus the designation if-diagonal is appro- 
priate. Another kind of identity is 

Pm{^(x)%^k'(x')} = 

unless K = K' and M = N + b. (57) 

where b is a vector in the same space as M and N, and 
whose only non-zero component is 6(K) = 1, so that 
M and N are identical apart from the change M(K) = 
N(K) + 1. 

There are many other relations like these, but all are 
essentially of the same kind, that the if-diagonal prop- 
erty is preserved only by matching a creation field oper- 
ator with a destruction field operator with the same K, 
though not necessarily the same x, either on the same or 
the other side of the density operator. If the matching is 
done on the same side, the configuration N is preserved, 
if the matching is done on opposite sides the configura- 
tion changes. 
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C. Derivation of the quantum kinetic master 
equation 



The equation of motion for the density operator is 



P = - % [H a 



H b + H e ,p]--[Hj,p] 



= C a p + C b p + C c p + C 2 p. 
Wc define the Laplace transform 



P(s) = 
and correspondingly 

Cn(s) = Vnp(s) 

W(s) = Qp(s) EE 

so that 

sv N (s) - i>n(0) 



t p(t)dA 



N 



(58) 
(59) 

(60) 

(61) 
(62) 



(C b + C c )v N (s) +Vtss£ 2 < ^«m(s) 

I M 



w(s)\ (63) 



and 

sw(s) — w(0) 



= {C a + C b + C c } w(s) + QC 2 J «m(s) 

I M 



+ w(s) } . (64) 



We assume w(0) = 0, that is the state that is already K- 
diagonal. This assumption cannot be made without some 
justification, which we shall postpone to Sect. IV D 3. 



Obviously we can choose any initial condition we wish; 
the justification required is that when we solve the equa- 
tions of motion it will be possible to show that w(t) 
rapidly becomes negligible, and negligible at all times, 
so that any time can be chosen as an initial time. 
From this we can then readily derive 



svn(s) - un(0) = (Cb + Cc)vn(s) + V-nC 2 vm(s) 

m 

+V N C 2 [s -C a -C b -C c - QC 2 ]- 1 QC 2 2 v M (s). (65) 



M 



This equation is exact; the physical content comes into 
the assumptions concerning what parts can be considered 
to be negligible. We invert the Laplace transform in (65), 
and consider the parts separately. 



1. Streaming and Quantum Terms 



These arc the names we use for C b and C c as defined 



in (J5J), (|41b|, |41c 

approximations 



). They keep the same form with no 
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2. The Forward Scattering Terms 



The term V(N)C2 X) ^m(s) can be simplified by not- 

M 

ing that %(s) is if-diagonal, and that 'P(N) projects 
onto a i^-diagonal density operator. Since C 2 is a com- 
mutator; i.e., 

£ 2 5> M (s) = -^X^' UM ( s )] (66) 

M M 

this can only yield a if-diagonal term from those parts 
of Hi which do not change the K-distribution; that is, 
the terms in Hi given by Hp, which we define as 

Hf = \ E ^(Ki,K 2 ,K!,K 2 ) 

Ki#K 2 

+ \ E tf(K 1 ,K a ,K 2 ,K 1 ) 

Ki#K 2 

+^C/(K 1 ,K 1 ,K 1 ,K 1 ). (67) 

This thus gives only transitions between the same K dis- 
tribution, and this can be regarded as forward scattering. 
These forward scattering terms give rise to the so-called 
mean field effects, that is the average effect on the motion 
of one particle of the interaction with all other particles. 



3. The Collision Terms 



We now introduce a Born approximation, which amounts 



to the neglect of Cb + £ c + Q£ 2 m the [ ] 1 term in (65) 
leaving a contribution to un(£) given by 



E 7 ^ 2 f 

M ''0 



dr exp{C a T}Q£ 2 v M (t - r). (68) 



Notice however that QC2 contains terms which depend 
on ipQ, ipQ, and when there is significant condensation 
these terms can become very large. Thus this approx- 
imation is valid only in the case of weak condensation. 
This means that we approximate the free evolution oper- 
ator by discretizing the eigenvalues to those correspond- 
ing to the centers of the momentum bands. This leads to 
an approximate measure of the validity of the neglect of 
the terms QC2 in comparison to the term retained, C a , 
which is given by requiring that the modification by the 
presence of any condensate to the excitation spectrum 
should be negligible for energy greater than h 2 A 2 /2m. 
Using the Bogoliubov theory, this leads to the conclusion 
that we must have the condition on the cell size 



where p is the density of particles. We now want to 
make a Markov approximation, which involves two steps. 
Firstly we neglect the r dependence in the last term, and 
secondly we let the upper limit in the integral become 
infinite. For this approximation to be valid there must be 
a smooth distribution over the available energy states, so 
that when we sum over all M, the resultant range of r in 
which the integrand is non-zero is very much smaller then 
the characteristic time over which vwi{t) evolves. This is 
a requirement on the kind of vyi{t) being considered. 

Now define the vector in the space of M, N, called 
e(1234) which can be written 

e(1234) Kl = e(1234) K2 = -e(1234) K3 = -e(1234) K4 = 1. 

(70) 

and all other components are zero, and let 

h 2 K 2 



2m 



We also introduce the notation 



AciJ e == OJ4. -\- UJ3 — LO\ — LO*2 , 



(71) 



(72) 



Before we examine the range of validity of them 
Markov conditions, we can see that using them ( |68[ ) now 
becomes 



1234 
l'2'Z' A' 



C/(l'2'3'4'), 
dr exp(i Acj e r)Q[[/(4321), ^ w m(*)] 



(73) 



Since the state umW is if-diagonal, and we project onto 
if-diagonal states, we must have 

1 = l',2 = 2' or 1 = 2',2 = 1' 
and 

3 = 3', 4 = 4' or 3 = 4', 4 = 3' 

but since the operator E/(1234) is symmetric in (12) and 
(34) this simply gives 4 times the result obtained from 
setting (1234) = (1'2'3'4'). 



Notice that no "forward scattering terms" arise in ( 73 ) . 
These have already been explicitly separated from the 
scattering terms and included in the term Hp defined in 
(|67j), and are eliminated from ( |73| ) by the Q projection 
operator. 

We can now carry out the time integral so that 
(fnl) becomes 



lim > 



h (e + iAui e ) 



(74) 
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Here we have used t/(1234) t = £7(4321). We now use 
1 1 



P- - inSiz) 
z + ie z w 

and rearrange, so that (|74|) becomes 



(75) 



e 

C/tC/e}. (76) 

e 

As is usual, we find a level shift term, which is a com- 
mutator, and a purely dissipative term. Adding together 
([76]), (66), (|5"9|), we get the full Quantum Kinetic Master 
Equation (QKME) 



i 

h 

i 

n 

i 

2h 



E / d3x (^) -JkW^nW 

K J 



ft 2 v 2 

2m 



^K(x),UN(t) 



53 (C/(K 1 ,K 2 ,K 1 ,K 2 ) + f/(K 1 ,K 2 ,K 2 ,K 1 )) +53[/(K 1 ,K 1 ,K 1 ,K 1 ),« N (t) 



Ki#K 2 



e ^ ' 

-pX]5(Aw(e)){2Cr(e)v N _e(t)l/ t (e) - U\e)U(e)v N (t) - v N (t)U^(e)U(e)}. 



(77a) 
(77b) 

(77c) 

(77d) 
(77e) 



D. Approximations 

1. Weak condensation condition 



Approximations only occur in the derivation of the 
collision and level shift terms. The major approxima- 
tion is the neglect of Cb + C c + Q£ 2 compared to C a . 
We choose A so that eigenvalues of C a (apart from zero 
eigenvalues) are substantially larger than typical mea- 
sures of the size of Q£ 2 on the particular states involved. 
Since the projectors do not affect the K = com- 
position of the states, the operator Q£ 2 includes the 
term (u/2) J (PxipQipQipoipo, which can become very large 
when there is significant condensation, hence this approx- 
imation will only be valid in the case of weak condensa- 
tion, the condition for which is given by ([59]); that is 



7r irh 
c A y/2mpu 



(78) 



We will deal with the issues of strong condensation in a 
subsequent paper. 
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2. The Markov Approximation 



The Markov approximation is really a kind of pertur- 
bative result. The change from the integral J Q dr — > 
Jq 00 dr requires the assumption that the terms that turn 
up have a broad spectrum of frequencies which are almost 
continuous. The frequencies which do turn up are of the 
form uj e , and the range of these is of the order kT/h, while 
a typical separation between the frequencies correspond- 
ing to adjacent transitions with a typical magnitude of 
momentum \j2mkT is of order of magnitude Ay // 2kT/m, 
which is the inverse of the time required for a particle to 
traverse a cell at a typical thermal velocity. 

The condition that we can regard the frequency spec- 
trum as almost continuous is that the separation between 
the frequencies is very much less than their range, that is, 
the density of states is sufficiently high. Quantitatively, 
this leads to 

l c > At = h/V2mkT. (79) 

If the density operators «n are such that the contri- 
butions from sums over the different e are all smoothly 
varying functions of the contributing frequencies, then 
the time dependence given by summing up over all tran- 
sitions in ( |73| ) will drop of over a characteristic time given 
by the inverse of the bandwidth, i.e, fi/kT, and this time 
must be very much less than the typical time scale of 
evolution of the the distribution, which is of the order of 
magnitude of the time between collisions. This leads to 
the requirement on the mean free path A m f p 

A mfp > A T . (80) 

Finally, the QKME itself gives a level broadening, 
which is implicitly small compared to the distance be- 
tween levels. The condition for this must be that there 
is little likelihood of a collision as a particle traverses a 
single cell, which means of course that 

A mfp > l c (81) 

All these conditions can be consistently satisfied. 

Coming back now to the quantum kinetic master equa- 
tion we see that the summations ^ e are over the discrete 
momentum ranges Ki,K2,K3,K4, and the energy con- 
servation delta function is sharp. This sharpness is really 
an artifact which arises from the extension of the upper 
limit of the time integral to infinity. This integral can 
be cut off at any time very much larger than the charac- 
teristic decay time of the kernel, Ti/kT . This means that 
the delta function can be taken as being broadened cor- 
respondingly, and it is easy to see that this broadening 
can encompass a significant number of momentum bands. 
Nevertheless, this broadened delta function can be writ- 
ten as an exact delta function if the remainder of the 
integrand is smooth, and the sums written as integrals 
using the continuum approximation ^2 ^ /^ 3 K- 

K 
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However, we have already assumed this smoothness in 
establishing the validity of the Markov approximation so 
this procedure is consistent with the Markov approxima- 
tion. 



3. The neglect of w(0) 

We will assume w(0) = 0, and then make an estimate 
of the size of w(t) and find under what conditions that it 



is negligible for all future times. From (64) it follows, on 
inverting Laplace transforms, that 

w(t)= f exp{(£ +£ 6 + £ c + Q£ 2 )(t-*')} 
Jo 

5>mW- (82) 



X 

M 



We can now approximate this by keeping C a , as in the 
derivation of the QKME, and replacing the effect of the 
all the other terms, including the interactions, by a phc- 
nomenological decay term, so that we write approxi- 
mately 



r e {(£.-7)N0}v^ ( i' )( i t ' 

JO M 



«;(*)« / el^-W-* n ^ VM (t')dt' (83) 

M 

= - { E E e (4A ^ 7)( ^ V(e), v m {t')]dt'. (84) 

e M 

As in the QKME derivation, we make the replacement 
Wjvr(i') «mW, and now we can do the integration over 
t' to get approximately 

w ^ ra ~\ E E ^-Me),VM (*)]. (85) 

e M ' 

For any given M and e we get a specific non if-diagonal 
term. The largest of these will occur when Aui e = 0. To 



estimate their size, we go to the definition (47) of U(e). 
Substituting the expansion on wavelet functions into this 
equation, and assuming the approximation it(x) = it<5(x), 
we find that the largest terms that occur are of the form 

J d 3 x/, v Kl (x, rJe^UK, (x, r))e iK -' x 

x^ 3 (x,r)) e - iK - x ^ 4 (x,r))e- iK4 - x 

xot(Ki, r)at (K 2 , r)a(K 3 , r)a(K 4 , r) (86) 

in which all wavelet functions refer to the same cell. Us- 
ing the definitions of the wavelet functions, the coefficient 
of the operator part in this equation can be estimated to 
be of order of magnitude The effect of the creation 
and destruction operators is of order of magnitude 1, so 
we can now estimate the size of the term in w(t) com- 
pared to those in Vjy[(i) to be of order of magnitude given 
by 
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, , ah 
w(t) « —j 



X 



E 

e.M,i 



7^ 3 

[at(Ki, r)at(K 2 , r)a(K 3 , r)a(K 4) r), v M (t)}. (87) 



The size of the initial coefficient can be made clearest by 
taking 7 to be given by VT/A m f p , in terms of which the 
coefficient can be simplified, leading to the condition for 
the validity of the neglect of non-if-diagonal terms: 

aAmfpAr ^ 



l 3 c 

This condition can be satisfied simultaneously with the 
two other conditions ( £9 ) and (|8(]). For example, for 
sodium at T = 2/iK, the scattering length is a = 4.9nm, 
leading to 

A m fp = 0.42m (89) 
A T = 4.8 x f(T 7 m (90) 
a = 4.9 x 10~ 9 m (91) 
(aA mfp A T ) 1/3 » l(T 5 m (92) 

Thus if l c is chosen in the range somewhat greater than 
10 _5 m this kind of treatment will be valid. 

The condition ( |69| ) is best viewed as a condition on the 
condensate density p; insertion of the sodium data into 
it gives the requirement that p <C 10 18 m~ 3 , which can be 
satisfied, but will not always be satisfied in current ex- 
periments. In fact this density corresponds to about one 
atom per cube of volume Ay using the above data, which 
indicates that for an accurate treatment of condensation, 
the weak condensation assumption will be invalid in this 
case. 

It should be borne in mind that all of the above has 
been in the absence of a trapping potential — the question 
will be revisited in our forthcoming paper dealing with 
strong condensation. 



4- Born Approximation 

The derivation of the collision terms relies on the Born 
approximation, and the collision rates appearing in the 
QKME ( |6a| ) are proportional to the Born scattering cross 
section. This implies that the potential should be weak, 
i.e. a small perturbation on the energy scale associated 
with the momentum coarse graining. In derivations of the 
quantum Boltzmann equation it is shown J5Q] that in a 
binary collision approximation, which assumes a low den- 
sity but not necessarily a weak potential (pa 3 < 1), the 
Born amplitude should be replaced by the two particle T- 
matrix element (proportional to the scattering amplitude 
a in low energy s-wave scattering). Essentially the same 
situation should pertain in this case and this means that 
we replace the exact interatomic potential with one for 
which the Born approximation is valid, and whose scat- 
tering length is the same as that of the exact potential. 
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However, a proof of this assertion in our formulation is 
at present lacking. 



5. Level shifts 

What are normally called "collisional level shifts" ap- 
pear in our formalism in two guises. There are the for- 
ward scattering terms which generate Hp , as defined in 
(|67j), and there are the higher order terms which from 



the derivation of the collisional terms in (77d). In the 
weak condensation situation which we want to treat in 
this paper, we expect the density to be so low that these 
can be ignored, but in our forthcoming work it will not 
be possible to do so without more careful justification. 



E. Stationary Solution 

Let us assume the level shift term is small — negligible 
in fact. Then it is clear, since none of the other terms can 
actually change the eigenvalue of H a , (the total coarse 
grained kinetic energy) or the total number of particles, 
that any function of 

^ = E^ a K(rW(r) (93) 

r,K 

and 

^ = E fl K(rW(r) (94) 

r,K 

will provide a stationary solution — clearly microcanoni- 
cal, canonical, and grand canonical versions will exist. 

Thus the result is the statistical mechanical result — 
but coarse grained in momentum space — the size A of 
the momentum coarse graining to be such that Hj is 
negligible compared with it on the states of interest. 



1. Correlation Functions in the Stationary Solution 

Since the stationary solution can be written as a func- 
tion of H a and N, the correlation functions of the wavelet 
creation and destruction operators can be written as 

(<4 c (r)o K '(r / )> =N K 6 K>K ,6 r>r ,. (95) 

Assuming the grand canonical form 

p,(T,n) = exp{-(H a - fiN)/kT}, (96) 

which is of quantum Gaussian form, we can easily see 
that 

(97) 



exp 



( flUJK 



\ kT 



- 1 



24 



The field operator correlation function is then 

(V^(x)^K'(x')} = <W.9(x-x')iVK(7» (98) 
and similarly 

(^k(x)^(x')> =«k,k'5(x-x')[1+JVk(T,m)]. (99) 

V. SIMPLEST APPLICATIONS OF THE 
QUANTUM KINETIC MASTER EQUATION 

The situation in which something like the "molecu- 
lar chaos" assumption of Boltzmann is valid is always of 
great interest, and it guides the intuition of most physi- 
cists in kinetic theory. We will consider in this section 
two simplified versions of the QKME which are valid in 
this situation, and one equation which permits the exis- 
tence of coherent effects, which cannot arise in a molec- 
ular chaos regime. 

The Uehling-Uhlenbeck equation, sometimes called the 
quantum Boltzmann equation, is a simple modification of 
the Boltzmann equation in which the phase space effects 
arising from quantum statistics are included. This equa- 
tion will be derived by assuming the distribution is not 
very different from equilibrium, and by factorizing higher 
order correlation functions. 

The situation in which we eliminate all dependence on 
position is obviously an extreme simplification, and cor- 
responds to the system being composed of only one spa- 
tial cell. However, in this situation the QKME reduces 
to a stochastic master equation of great simplicity and 
intuitive appeal. This equation does not appear ever to 
have been treated in the literature, so we have given it 
the name "quantum Boltzmann master equation" . 

Finally, a very simplified quantum mechanical master 
equation for the K = band can be derived by assum- 
ing all other bands are thermalized. This equation will 
be called the "condensate Master equation" , and we are 
able to give an appealing description of the initiation of 
the condensation process, and of the growth of the con- 
densate as long as the weak condensation condition is 
valid. 



A. The quantum Boltzmann master equation 



A very simple stochastic master equation can be de- 
rived by assuming that I — ir/ A is equal to the size of the 
whole system, so that r takes on only one value. In this 
case we find that the K-bands a re al l one- dimensional, so 
that all the commutator terms, (77a- 77d ) are zero. Since 
in this case there is really no difference between n and 
N, we can write 



v n (t) = |n)(n|w n (t) 



(100) 
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so that w n (t) is the occupation probability of the state n. 
We then derive what we will call the quantum Boltzmann 
master equation 



n = X! S ^ LJi + ^3 - ^2 - ^l)|C/l234| 2 

1234 

x|nin 2 (n 3 + l)(n 4 + l)[w a - w n+e ] 

+ {n 1 + l)(n 2 + l)n 3 7i 4 [u; n - w n -e] j (101) 



For clarity of argument we will also write (100) in the 
form 



w(n) = {*i234( n - e ) w ( n - e ) - *i234( n >( n )j 

1234 

+ {wC" + e)w(n + e) - t U34 {n)w{n)} (102) 

1234 

where i^^l 11 ) = unless uj± + lu2 = cj 3 + w 4 , and if this 
is satisfied 

i i234( n ) = 7i234(ni + l)(n 2 + l)n 3 n 4 

*r234( n ) = 7i234nin 2 (n 3 + 1)(«4 + 1) (103) 



and 



71234 = 72 ■ (104) 



T/ie Boltzmann master equation 



If we drop the terms 1 + in the quantum Boltzmann 
master equation we get an equation called the Boltzmann 
Master Equation, which has been previously considered 
by Van Kampen (m) and Gardiner |55fl based on stochas- 
tic arguments. The equation we have given includes no 
streaming terms, since we have only one spatial cell, and 
thus our treatment assumes that the energy levels are dis- 
cretely spaced. This means that our normalization vol- 
ume is so small that the broadening of the levels which 
arises from collisions is significantly less than the spacing 
between levels. 



2. Stationary solutions 



Apart from the ambiguities caused b y con served quan- 
tities, the solutions of equations like (101) arc unique. 
There is a consistent detailed balance stationary solu- 
tion, in which 



*i234( n - e Wn - e) = t 1234 (n)w s (n) 



(105) 



and from this we see, using the explicit form ( |103| ) of t^, 
that 
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w s (n-e) _ t 1234 (n) 



(106) 



so that for all n able to be connected by a collision 1+2 ^ 
3 + 4, 

w s (n) = constant. (107) 

Thus w s (n) is a function only of conserved quantities, 
total kinetic energy, momentum and the total number 
of particles. We can thus choose the grand canonical 
solution given by 

w s (n) ex exp f >t— J (108) 

where E is the total kinetic energy, N the total number 
of particles, and P the total momentum of the particles, 
while u is the mean velocity of the system. The station- 
ary distribution is in this case factorizable 

to a (n)=JJ(e ) (109) 



3. Factorized equation for the mean occupation numbers 
We can straightforwardly derive the equation for (n a ) = 

(h a ) = 4^7 a23 4| - (nan 2 (n 3 + l)(n 4 + 1)) 

234 

+ ((n a + l)(n 2 + l)n 3 n 4 )} (110) 



If we assume that the averages inside this equation can 
be factorized, which should be valid if we are not too far 
from equilibrium, and for compactness set (rij) — ► rii, we 
get (without the streaming terms) the Uehling-Uhlenbeck 
equation (42); 

h a = 4^7 a 234| - n a n 2 (n 3 + l)(n 4 + 1) 

234 

+(n a + l)(n 2 + l)n 3 n 4 } (111) 

The quantum Boltzmann master equation is a very sim- 
ple equation whose simulation is entirely feasible, and 
results of such simulations will be presented elsewhere. 
However, as shown by our derivation, it is expected that 
it will not give an accurate description of situations in 
which there is significant condensation. Nevertheless it 
could well give valuable insights into the region just below 
the threshold of condensation, which will be investigated 
in QK1TL 



27 



B. The Uehling-Uhlenbeck equation 



We will now show that we can deri ve from the quan- 
tum kinetic master equation ( 77a - 77e ) the more conven- 
tional kinetic equation known as the Uehling-Uhlenbeck 
equa tion |42j |. This equation is essentially the same as 
( 111 ), but with streaming terms added as well, which in 
our derivation will arise naturally out of the streaming 
terms of the QKME. For this we consider the one particle 
distribution function to be defined as 



/k(x) = (J) 3 THpV4c(x)Vk(x)} 

= (x) 3 ^ Tr{l,N ^( x )^(x)}. 

N 

Notice that the definitions ( |3l| , |34| ) imply that 
tMx) = 5>-* k °Vk(x) 

K 



(112) 
(113) 

(114) 



so that the usual probability density function is, because 
the wn are if-diagonal, 



-) Tr{pV t (x)V(x)}=^/ K (x 



(115) 



so that /k (x) can be viewed as a phase space density per 
phase space volume h 3 . 



1. Streaming terms 

We consider the various terms in the quan tum kinetic 
master equation. The term arising from ( [77a| ) is obtained 
by considering 



V>k(x), / d 3 y^K(y) 
which can be written 

V;k(x),| d 3 y^(y) (™) VVk(v) 
= J d 3 y 5 (x-y) f-^W^ K (y 
and now integrating by parts, 

dill) = -^V x (J d 3 yg{ X - y)^ K (y 



(116) 

(117) 
(118) 

(119) 



but the integral is simply the projection of "0k (y) on to 
the K subspacc; thus 



O = VVk(x). 



(120) 



2S 



Using this kind of technique, it is straightf orwa r d to 
show that the contributions to /k(x) from (77a, f77ej ) 
become 



/k(x)U 



KK 



V x /(K,x) 



/ K (x)| fc =(J) 3 V-j K (x) 



(121) 



(122) 



where 



ih 



Jk(x) = — (v£(x)VVk(x) - [V^(x)]Vk(x)) . (123) 
2m \ / 



The term ( 121 ) is the usual Boltzmann streaming term. 
The term (122) is a purely quantum mechanical term 
arising from the spreading of the wavepacket. 

In order to get some idea of the significance of this and 
some terms to come, we make a local equilibrium ansatz, 
in which we write 



(^(x)^(x')>=3(x-x')/, 



K 



X + x' 



?K,K' 



(124) 



This is a natural generalization of (|98|) to a slightly non 
equilibrium situation, and since g(0) = (7r/A) 3 , it is con- 



sistent with (112). Using this form we find quite straight- 
forwardly that Jk(x) = 0. The extent to which this term 
does not vanish must be related to the extent to which 
the mean momentum in the band (K — A, K + A) is 
different from TiK, and this difference is of order of mag- 
nitude A, which will only be appreciable for very small 
K, and in particular for K = 0. 



2. Forward scattering terms 
If we work out the mean value of the commutator, we 



find for the term arising from (77c) 

-{2 J d 3 x'X](V^(x')^(x')V'K'(x')^K(x)). 9 (x-x'). 
-2 ( d 3 x'(^(x)V4(x')^K'(x')^ K (x')) 5 (x-x')}. 



I 

Ti 



(125) 

If g(x — x') was a perfect delta function, this term would 
vanish, thus the degree to which it does not vanish de- 
pends on the smoothness of the averages as functions of 
x, x'. 



Using the local equilibrium ansatz (124), and a Gaus- 
sian factorization assumption on the four point correla- 
tions, it is easy to show that this term vanishes, and is 
therefore presumably very small when the state is close 
to equilibrium. Thus we can say /k(x)| c ~ 0. 
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3. Collision terms 



In evaluating the contribution from (77e) we will get a 
number of terms, of which a typical one is proportional 
to 



J2 I d 3 y d 3 y' 5 (x - /) (lb (x)^ 2 (y>K 3 (y')V>K 4 (y')^ 4 (y)^ 3 (y) ^ 2 (y)< (*)> 

K 2 ,K 3 ,K 4 J 

x^ K + ^k 2 - wk 3 - WK4 ) e *( K + K 2 - K ^^)-(y-y'). (126) 



If we factorize as in local equilibrium, and use (124), we 
get terms like 



y' + y 



E / rf 3 ydV.9(x-y')g(x-y)[.g(y'-y)] 3 /K(^)/K 2 

K 2 ,K 3 ,K 4 ^ ' 

x[/k 3 (^iA+W^ f^") +l]^K + ^K 2 -^K3-^K 4 )eW- K 3-K 4 ) (y-y'). (127) 



Assuming the functions are relatively smooth, so that 
we can set (because of the peaked nature of g(x — y'), 
ff(x-y) and.g(y'-y) 



^«^x (128) 



we get 



E /K (x)/k 2 (x) [/k 3 (x) + 1] [/k 4 (x) + 1]S(UJ K + UJ K2 - WK 3 - WK 4 ) 
K 2 ,K 3 ,K 4 

x|rf 3 y| rf 3 y'.9(y')5(y)[3(y'-y)] 3 e i(K+K2 - K3 - KlHy - y ' ) . (129) 

The last factor is in fact able to be evaluated exactly, 
and the result is a momentum conservation function 
M A (K + K 2 - K 3 - K 4 ), where 

M A (Q)=(-) Hq6 Qi , + l5 Qi , A + l5 Qi ,- A }. (130) 
^ n ' i=i 

Thus "momentum" K is not exactly conserved, since the 
quantity K in i/>k(x) represents only the midpoint of the 
range of possible momenta; it is possible for three mo- 
menta, each from such a range, to add up to a value 
either within the range, or in either of the ranges above 
or below. The exact proportions of each are shown by 
the calculation to be | : g : ^. However, assuming that 
/k(x) are smooth functions of K on the scale of A, we 
can ignore the slight spreading, and make the transfer 
from a summation to an integration: This yields 



E (-) m a(K + K 2 - K 3 - K 4 ) -> — ^ f d 3 K 2 f d 3 K 3 f d 3 K 4 f 6(K + K 2 - K 3 - K 4 ) 

K 2 ,K 3 ,K 4 \ 7r ^ ( J J J J 

(131) 

so that this form eventually becomes 
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-^2 J d 3 K 2 J d 3 K 3 J d 3 K 4 6(K + K 2 - K 3 - K 4 )<5(^ + lo 2 - u 3 - w 4 )/k(x)/k 2 (x)[/ Ks (x) + 1][/k 4 (x) + 1] 

(132) 

Here we have made the approximation 

u(x - x') = u5(x - x') (133) 
Taking all terms into account we eventually get a term 

/(K, x)| e = y d 3 K 2 J d 3 K 3 J d 3 K 4 6(K + K 2 - K 3 - K 4 )<^ + u 2 - u; 3 - « 4 ) 

x{/ k (x)/k 2 (x)[/k 3 (x) + 1][/k 4 (x) + 1] - [/ K (x) + 1][/k 2 (x) + 1]/k 3 (x)/k 4 (x)} (134) 

which is the collision term of the equation known as the 
Uehling-Uhlenbeck equation. 

4- Summary 

In the situation that we have an approximate local 
equilibrium, and that multifield correlation functions can 
be factorized as if the density operator is Gaussian, we de- 
rive the Uehling-Uhlenbeck equation as an approximate 
kinetic equation. 

^(K, x) « ^ K m Vx /(K, x) + j d 3 K 2 J d 3 K 3 J d 3 K 4 (5(K + K 2 - K 3 - K^lu + u; 2 -lu 3 - lu 4 ) 
x{/k(x)/k 2 (x)[/k 3 (x) + 1][/k 4 (x) + 1] - [/ K (x) + 1][/k 2 (x) + 1]/ Ks (x)/k 4 (x)} 



(135) 



In situations of interest for Bose Condensation it would 
no longer be valid to make the factorizations necessary 
for the collision term as above, nor factorization in the 



local equilibrium assumption (124) 



C. Equations for the condensate density operator 

Let us now assume that the scattering of non conden- 
sate modes is strong, and that we are only interested in 
the behavior of the condensate. We can then assume that 
the total density operator can be factored into a thermal 
non-condensate part, and a condensate density operator 
Po(t) as a first approximation at least. Thus 

p(t) = ^v N (t)->p B (ji,T)®p (t) (136) 

N 

This situation is not that employed in current experi- 
ments, where the condensate grows by taking atoms out 
of the bath of warmer atoms, but is given for illustrative 
purposes only — it describes the situation in which the 
condensate mode is put in contact with a bath composed 
of all the atoms in the other modes, which are held at a 
given temperature and chemical potential. More realistic 
treatments of condensate growth are left to a later paper. 
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We trace over the non-condensate modes and the as- 
sumption that the non-condensate is thermalized will 
mean that all terms not involving condensate operators 
in the QKME will be vanish. The only non- vanishing dis- 
sipative terms will involve matched ipQ and ipo operators, 
so that we arrive at the condensate master equation: 



Po(t) = 



d 3 x^(x) 



h 2 v 2 



1. 



Vo(x) + -C/(0000),poW 
2m I 2 



(ug(0) «k) 



d 3 xipl(x)ip (x),p 



+ I d 3 xj d 3 x'G(-)(x-x',T, M ){2^o(x)p V'S(x')-/OoV'S(x')V'o(x)-^(x , )Vo(x)po} 

d 3 x J d 3 x'GW(x-x' ) T, / i){2^(x) /9o Vo(x / )-po^o(x / )4W-^o(x / )^(x)po}, (137) 
in which the quantities > are given by 

G«(x - x', T, p) = Tr B jg ^Zt(l, 2, 3, x)Z(l, 2, 3, x')ps(T, p) j 
G(-) (x - x', T, p) = Tr B J ^2 Z (h 2, 3, x)^t(l, 2, 3, x')pb(T, p) 1 . 

I 123 ^ J 

(138) 

with the definition 

Z(K 1; K 2 ,K 3 ,x) = ue - 4 ( Kl+K2 - K3 )- x ^ Kl (x)^K 2 (x)^ 3 (x) (139) 
Using thermal averages we find 

GW(x-x',T,p) = ^{^[ 9 ( x -x / )] 3 e i ( Kl+K2 - K ^-( x - x ')rlK 1 nK 2 (nK3 + +w 2 -^ 3 )} 

" 123 
2 

G(-)(x-x',T,p) = ^{^[ 5 ( x -x')] 3 e l ( Kl+K2 - K3 H— ')(n Kl + l)(n Ka + l)n K3 <5(wi + w 2 - w 3 )} (140) 

" 123 



Suppose we now set, as in (|97[ ) 
then it follows that 



(K ^0) (141) 



^'^(1 + ^3) = e M/feT (n Kl + l)(n K2 + l)n K3 , (142) 
so that 

G (+) (x - x', T, p) - e^ fcT G ( - } (x — x',T, p). (143) 



Quantum stochastic differential equation form 

We can write a Quantum stochastic differential equa- 
tion (QSDE) equivalent to the master equation (137) by 



using the methods in Chap. 5 of Jl0|]. This equation takes 
the form 
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ihdipo(y,t) = | - ^VVo(y>*) - \ ug(Q) n K j tp (y,t)+u J d 3 x.g(x - y)V>o( x > t)tp (x, t)ip (x, t) 

+ J d 3 x J d 3 x' 3 (x'-y){GW(x-x',T, M )-G(-)(x-x',T, M )}^ (x,i)| 

+ y d 3 xg(x-y)dW{x 7 t) (144) 



in which dVK(x, t) is a quantum white noise increment 
which satisfies the conditions 

dW(x,i)dWt(x',t) = 2G ( - ) (x-x',7» (145a) 

dM^ t (x,t)dVK(x',i) = 2G (+) (x-x',T jA t) (145b) 

dW(x, i)dW(x / , t) = (145c) 

dW t (x,i)dW t (x',t) = (145d) 



2. Gain and ioss 



In order to get some idea of the predictions of the mas- 
ter equation (137) we can write equations for the averages 



(^o(y)) = ^(y), 



(146) 
(147) 



These become 



aft(y) 
<9i 



gp(y) 
at 



— V 2 0(y) - u fl (0) J] "K 0(y) + u / d 3 x 3 (x - y)(v4(x)Vo(x)V>o(x)> 

+ |d 3 x| d 3 x'. 9 (x'-y){G( + )(x-x',T, A1 )-G(-)(x-x',T, / ,)} ( /)(x) 
-V-j(y)+ J d 3 x{G l+) (y-x,T,/i)-G ( ->(y-x,T, /1 )} 
x(^(x)^ (y)+^(y)Vo(x)) + 2 / rf 3 x G(+)(x,T, M ) 5 (x). 



Here j(y) is the probability current. The equation (148) 
for 4>{y) is a rather familiar non-linear Schrodinger equa- 
tion form, but contains the second line as well, which is 



an explicit gain- loss term. We can use the relation (142) 
to write this in the form 



d 3 x / d 3 x'g(x'-y)( e ^ T -l)G(+)(x-x',T, M )0(y) (149) 



and clearly when fi = 0, this vanishes. 

If we consider spatially homogeneous solutions, then 
we can write </>(x) — > 0, and we find 

^ = -^(ug^^flK^^- J d 3 x 5 (x-y)(v4(x)V>o(x)^o(x)) 



(148) 



+4>{e^ /kT - l)-4-r / d 3 Ki / rf 3 K 2 n Kl n K2 (l + "k 1 +k 2 )'5(wk 1 +w K , -w Ki +k 2 )- 
2n it J J 

(150) 
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3. Stationary solutions. 



If fi < 0, as must be the case above the condensation 
temperature, the last part clearly represents a loss term. 
There will be no growth of a non zero <f> value. If we 
assume a Gaussian solution, we can write 

(4(y)My)My)) W4(y)My)) (i5i) 

and it is clear that the only stationary solution will be 
= 0. 

Sim ilarl y, we can get an approximate closed equation 
from fll4g| ) by assuming 

(^(x)Vo(y)) = £ ^5p* ( 152 ) 

which assumes the equilibrium shape ( p9| ) of the correla- 
tion function. This leads to 

g=2 [d 3 xG+(y-x,T >t j,)g(x-y){(e^ kT - l)p/g(0) + 1}. (153) 
We then find the stationary solution 

P = i _ ^J/kT = I _ e n/kT > ( 154 ) 

which gives the total mean number per volume as / 3 as 
1/(1— e^/ kT ), the usual statistical mechanical result when 
/i < 0. Whe n p, = there is no longer a stationary state 



of Eq..(148); p will grow indefinitely. However there is 
still no gain, so no coherent phase appears. Of course in 
practice fx is not exactly zero, so indefinite growth will 
not occur. 

Thus we see that the ideal Bose condensate result arises 
purely out of noise; that no coherent phase appears at 
all. This result is only valid in the weak condensation 
limit; we will see in a subsequent paper that higher order 
terms can lead to an effective gain, and that a non-zero 
stationary value of 4> can arise. 



4 ■ Gain and a coherent phase arising from a modified 
distribution of energies in the non- condensate 

It is of interest to see whether we could produce a non- 
zero value of 4> by manipulating the nature of the energy 
distribution of the non-condensate. If we do assume that 
the non-condensing modes are in a non-equilibrium sta- 
tionary state, so that instead of pb{T,/i) we h ave an ar- 
bitrary density operator p ne , the gain in ( |148| ) becomes 
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u 2 



2h 2 



— J d 3 Kid 3 K 2 (5(a) Kl +WK 2 - ^K 1 +K 2 )|riK 1 nK 2 [l + n Kl +K 2 ] - [1 + "Ki ] [1 + ™K 2 ]n Kl +K 2 | • 

(155) 



We can get gain if (but not only if) every team in the {} is positive. Assuming that n(K) is only a function of cjk, 
this will happen when 

n(wi) n(w 2 ) l + n{uj 1 +uj2) >x 



1 + n(u>i) 1 + n(uj 2 ) n(uji+uj 2 ) 
or, defining 

F(uj 1 ) + F(lu 2 )> F(lo 1 +lu 2 ). (158) 
This is the condition that F(v) is a convex downwards function as in the solid line in the diagram. For example, 

* . — hu — Xuj 2 , . 

FH = — kf — (159) 



gives 



and the net gain is 



F(u> 1 )+F(u 2 )-F(u 1 +u 2 ) = ^p^>0 (160) 



u 2 f 
2h ir b J 



xcxp( — ) n Kl n K2 [l + n Kl +K 2 ]- (161) 



VI. CONCLUSIONS 



Quantum kinetic theory is both a genuine kinetic theory and a genuine quantum theory. The kinetic aspect arises 
from the decorrelation between different momentum bands, and this is essentially our version of the molecular chaos 
condition. However, because each K-band has a range of momenta within it, from which coherences are not eliminated, 
there is still the possibility of long wavelength quantum coherence. In practice such coherence is not expected to be 
of significance except in the Bose condensate, but the formulation we have developed enables a correct formulation of 
the method by which incoherence can be transferred into the condensate. 

We have formulated quantum kinetic theory in terms of the Quantum Kinetic Master Equation for bosonic atoms. 
It is a quantum stochastic equation for the kinetics of a dilute quantum degenerate bose gas, and it is valid in the 
regime of Bose condensation, as well as in regimes where no condensation takes place. 

The QKME is a genuine TV-atom equation, which is intermediate between the full description in terms of the 
complete density matrix and the familiar kinetic equations for single particle distribution functions (like the quantum 
Boltzmann equation). It contains as limiting cases both the Uehling-Uhlcnbcck (or quantum Boltzmann equation) 
and the Gross-Pitaevskii equation. 

The key assumption in deriving the QKME is a Markov approximation for the atomic collision terms. This is valid 
to the extent we are interested in a situation where in addition to the Bose condensate we have a thermal fraction of 
atoms which effectively acts like a heat bath with a (short) thermal correlation time H/kT, which is taken as zero in 
the Markov approximation. 

The present paper has developed the basic structure of the theory, has stated and justified the approximations 
which are needed, and has thus delineated the region of validity of the theory. Extensions to include a trapping 
potential and to account for a large fraction of condensed atoms giving rise to Bogoliubov type excitation spectrum 
are essentially technical, because the basic formulation carries through with little conceptual change, although the 
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equations which are appropriate in these cases are significantly changed. These aspects, which are essential for the 
application of the QKME method to realistic experiments, will be dealt with in QKII 56 1. 
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[1] M. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Science 269, 198 (1995). 
[2] C.C. Bradley, C.A. Sackett, J.J. Toilet, and R. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 

[3] K. B. Davis , M.-O. Mewes , M. R. Andrews , N. J. van Druten , D. S. Durfee , D. M. Kurn , W. Ketterle, Phys. Rev. Lett. 
75, 3969 (1995); 

[4] M.-O. Mewes , M. R. Andrews , N. J. van Druten , D. M. Kurn , D. S. Durfee , C. G. Townsend , W. Ketterle , Phys. Rev. 
Lett. 77, 416 (1996). 

[5] D.S. Jin, J. R. Ensher, M.R. Matthews, C. E. Wieman, E.A. Cornell, Phys. Rev. Lett. 77, 420 (1996). 

[6] M.-O. Mewes , M. R. Andrews , N. J. van Druten , D. M. Kurn , D. S. Durfee , C. G. Townsend , W. Ketterle, Phys. Rev. 
Lett. 77, 988 (1996). 

[7] D.W. Snoke, J. P. Wolfe, and A. Mysyrowiccz, Phys. Rev. B 41, 11171 (1990); J.-L. Lin and J. P. Wolfe, Phys. Rev. Lett. 
71, 1222 (1993). 

[8] Yu. M. Kagan, B. V. Svistunov and G.V. Shlyapnikov, Sov. Phys JETP 75, 387 (1992). 

[9] Yu. Kagan and B.V. Svistunov, JETP 78, 187 (1994). 
[10] Yu. Kagan, G. V. Shlyapnikov, J. T. M. Walraven, Phys. Rev. Lett. 76, 2670 (1996). 
[11] O. J. Luiten, M. W. Reynolds, J. T. M. Walraven, Phys. Rev. A 53, 381 (1996). 
[12] Yu. Kagan, E. L. Surkov, G.V. Shlyapnikov, Phys. Rev. A 54, R1753 (1996). 
[13] H.T.C. Stoof, Phys. Rev. Lett. 66, 3148 (1991). 
[14] H.T.C. Stoof, Phys. Rev. A 49, 3824 (1994). 

[15] H. T. C. Stoof, M. Houbiers, C. A. Sackett, R. G. Hulet, Phys. Rev. Lett. 76, 1 (1996). 

[16] E. Levich and V. Yakhot, Phys. Rev. B 15, 243 (1977); J. Phys. A 11, 2237 (1978). 

[17] D.W. Snoke and J.P. Wolfe, Phys. Rev. B 39, 4030 (1989). 

[18] D.V. Semikoz and I.I. Tkachev, Phys. Rev. Lett. 74, 3093 (1995). 

[19] V.V. Goldman, I.F. Silvera, and A.J. Leggett, Phys. Rev. B 24, 2870 (1981). 

[20] D.A. Huse and E.D. Siggia, J. Low. Temp. Phys. 46, 137 (1982). 

[21] R.V.E. Lovelace and T.J. Tommila, Phys. Rev. A 35, 3597 (1987). 

[22] P. Ruprecht, M.J. Holland, K. Burnett, and M. Edwards, Phys. Rev. A 51, 4704 (1995). 
[23] M. Edwards and K. Burnett., Phys. Rev. A 51, 1382 (1995) 
[24] R. Ballagh and K. Burnett, submitted to Phys. Rev. A 

[25] M. Edwards, R. J. Dodd, C. W. Clark, P. A. Ruprecht, K. Burnett, Phys. Rev. A 53, R1950 (1996). 

[26] Mark Edwards, P. A. Ruprecht, K. Burnett, R. J. Dodd, and Charles W. Clark, Phys. Rev. Lett., Phys. Rev. Lett. 77 , Aug 
(1996). 

[27] G. Baym and C.J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

[28] S. Stringari, Phys. Rev. Lett. 76, 1405 (1996). 

[29] F. Dalfovo, and S. Stringari, Phys. Rev. A 53, 2477 (1996). 

[30] S. Stringari, Phys. Rev. Lett 77, 2360 (1996). 

[31] LYou and M. Holland, Phys. Rev. A 53, 1 (1996). 

[32] M.J. Holland and J. Cooper, Phys. Rev. A, 53, R1954 (1996). 

[33] M. Holland, J. Williams, K. Coakley, and J. Cooper, Quantum Semiclass. Opt. 8, 571 (1996). 
[34] M. Lewenstein and L. You, Phys. Rev. A 53, 909 (1996). 
[35] L. You, W. Hoston, M. Lewenstein, and K. Huang, preprint. 
[36] R. Dum, Y. Castin and J. Dalibard, preprint. 
[37] B. Esry, preprint. 

[38] Victor M. Perez-Garcfa, H. Michinel, J.I. Cirac, M. Lewenstein, and P. Zoller, accepted in Phys. Rev. Lett. 
[39] A. Griffin, Phys. Rev. B 53, 9341 (1996). 

[40] C. W. Gardiner, Quantum Noise (Springer- Verlag, Berlin Heidelberg, 1991). 

[41] Bose-Einstein Condensation, edited by A. Griffin, D. W. Snoke, and S. Stringari (Cambridge University Press, Cambridge, 
1995). 

[42] L.W. Nordheim, Proc. Roy. Soc. A, 119, 689 (1928); F. Bloch, Z. Phys. 52, 555, (1928); (Initial derivations of the Uehling 
Uhlenbeck equation); S. Kikuchi and L. Nordheim, Z. Phys. 60, 652 (1930); (More rigorous derivation of the Uehling 



36 



Uhlenbeck equation); E. A. Uehling and G.E. Uhlenbeck, Phys. Rev. 43, 552 (1933) (Application of the Uehling Uhlenbeck 
equation to transport phenomena in quantum gases.) L. P. Kadanoff and G. Baym, Quantum statistical mechanics : Green's 
function methods in equilibrium, and nonequilibrium (Addison- Wesley, 1962). See also S| for a more recent derivation, and 
references to other contemporary work. 
[43] V.L. Ginzburg and L.P. Pitaevskii, Zh. Eksp. Teor.Fiz. 34, 1240 (1958) [Sov. Phys. JETP 7, 858 (1958)]; E.P. Gross, J. 
Math. Phys. 4, 195 (1963). 

[44] P.D. Drummond and S.J. Carter, J. Opt. Soc. Am B, 4, 1565 (1987); S.J. Carter and P.D. Drummond, Phys. Rev. Lett. 67 
(1991) 

[45] A. A. Abrikosov, L.P. Gorkov and I.E. Dzyaloshinski, Methods of Quantum Field Theory in Statistical Physics, Dover NY 
(1963) 

[46] E.M. Lifshitz and L.P. Pitaevskii, Statistical Physics Part 2, Landau and Lifshitz Course of Theoretical Physics Vol. 9 
(Pergamon Press, Oxford 1980) 

[47] J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 68, 580 (1992); C. W. Gardiner, A. S. Parkins, and P. Zoller, Phys. 

Rev. A 46, 4363 (1992); H.J.Carmichael, An Open Systems Approach to Quantum Optics, Lectures Notes in Physics ml8 

(Springer, Berlin, 1993). 
[48] D. Jaksch, C.W. Gardiner, and P. Zoller, unpublished. 

[49] A.J. Leggett, in Bose-Einstein Condensation, Eds. A. Griffin, D.W. Snoke and S. Stringari (Cambridge University Press, 

Cambridge 1995), p. 452. 
[50] see, for example, Ref. 

[51] S. N. Bose, Z. Phys. 26, 178 (1924), A. Einstein, Sitzungsber. Kgl. Preuss. Akad. Wiss. 1924, 261 (1924). 
[52] J. Ross and J.G. Kirkwood, J. Chem. Phys. 22, 1094 (1954) 

[53] C. K. Chui, An introduction to Wavelets, Academic Press (Boston 1992); Y. T. Chan Wavelet Basics, Kluwer Academic 
(Boston 1995) 

[54] N. G. van Kampen, Stochastic Processes in Physics and Chemistry, North Holland (Amsterdam 1981,1983) 
[55] C.W. Gardiner, Handbook of Stochastic Methods (Springer, Heidelberg, 1983, 1985, 1989). 
[56] C.W. Gardiner and P. Zoller, unpublished. 



37 



K " 




Kl 



K 2 



•••• 



N(Ki) 



N(K 2 ) 



K=0 



K 

(0 




K 2 




K 3 
K 4 

K 2 



x)u{\r-x\) 



Vg(x,r) 



